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The one-dimensional Kondo lattice model is investigated by using bosoniza- 
tion techniques and conformal field theory. In the half-filled band, the charge and 
spin gaps open for the anti-ferromagnetic Kondo coupling. Away from half-filling, 
the paramagnetic metallic state is characterized by the fixed point of Tomonaga- 
Luttinger liquid with the large Fermi surface in accordance with known results. It is 
suggested that as the electron density approaches half-filling, both of spin and charge 
susceptibilities may show a divergence property. 
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There is much current interest in the Kondo lattice model (KLM) , which is considered to 
be a basic model for heavy fermion systems. The model hamiltonian consists of conduction 
electrons coupled with a localized spin array via the Kondo exchange interaction. The 
competition between the Kondo effect and the RKKY interaction results in various phases 
such as magnetic phases, Kondo insulators, etc. As a first step to understand the KLM, the 
one-dimensional (ID) KLM has been studied extensively by renormalization group methods 
IHH], exact analytic methods ||, numerical diagonalizations [j§Jj]], Monte Carlo simulations 
fU?], and bosonization methods f|||, which have clarified basic properties of the model. In 
this paper we investigate the low-energy physics of the ID KLM by using bosonization and 
conformal field theory (CFT) techniques, and give complementary discussions to the results 
known so far. In particular, we point out that the marginally relevant spin interaction plays 
a key role in the model. 

We consider the ID KLM, 

H = -t^c\ a c i+lia + \ K ^S c ,i ■ S fti + X f ^2Sf ti ■ Sf ii+ x, (X K ,X f >0), (1) 

i,a i i 

where the interaction A/ for /-spins is introduced, which makes it easy to treat localized 
spins by bosonization methods. Here Cj )CT and c\ a , are the annihilation and creation operators 
of conduction electrons, and S c = c] Q ,<T Q!)( gCi )J g/2 with the Pauli matrix cr. Similarly, Sf = 
fia (T <x,pfi,/3/2 for localized /-electrons. In order to study the low-energy critical behavior, 
we use bosonization techniques combined with CFT. Applying a non-abelian bosonization, 
we first separate the charge and spin degrees of freedom of conduction electrons preserving 
SU(2) symmetry of the spin sector [H^-12|. This method enables us to observe whether the 
interaction is relevant for the gap formation. In a continuum limit, the electron operators are 
expressed in terms of the left-going (L) and right-going (R) operators, c a (x) = e lkpl ''CL a (x) + 
e~ tkFX CR a (x). Introducing the current operators for the charge and spin degrees of freedom, 
Jc,l = c^ La CLa, J c,l = c^ La ^(T a pCLi3, etc., we can represent the hamiltonian by these currents, 

H = H c + H s , 

H c = — J dx[J C)L (x)J C)L (x) + J c ,r{x)J c ,r{x)1 

2 



H s = J dx\J C)L (x) ■ J c ,l(x) + J C)R {x) ■ J c ,r{ x )} + )kk J d xS c(x) ■ Sf(x) 

+ / MJfA%) ■ JfA%) + JfA%) ■ JfAx)}, (2) 

where v c (v) and Vf are the velocities of the spin (charge) excitation of conduction electrons 
and localized /-electrons. J c ,l(r) an d J c (f),L(R) satisfy U(l) current algebra and level- 1 SU(2) 
current algebra [10- E|. From non-abelian bosonization formulas, c\ a CRp oc (g c )ap£ l "^^ c > 
etc., we can express the spin operator as 

S c = J C , L + J C , R + const.(e 2 ^tr(£ c oV^ C) + h.c), (3) 

where g c is the fundamental representation of SU(2) xSU(2), and is the bosonic phase 
field related to the charge degrees of freedom. A similar formula holds for the localized spin 
S f with the charge part exp (iv / 2vr0j' ) ) being replaced by its expectation value. 

Let us begin with the half-filled case in which the number of conduction electrons equals 
that of lattice sites. In this case, it has been deduced numerically that the excitation gaps 
open both in the charge and spin sectors characterizing the Kondo insulator Q. We shall 
discuss this analytically below (see also |J). The Kondo interaction in eq.(@) is rewritten in 
terms of boson fields, 

H^t = Ajf / dx(J C:L (x) + J c>R (x))(J f>L (x) + J fA x )) 

+ ^const. J dx[tr(^ c cr)tr(^ / o-)e iv/ ^^ C) + ti(g c a)ti(g\a)e iV ^^ C) + h.c.}, (4) 

where we have dropped the irrelevant oscillating terms. Before switching on the exchange 
interaction, the scaling dimensions of the fields exp(i\/27r0 c ), g C) and gf are all equal to 



1/2 fl3| , |14 |. Therefore the second term of eq.@ turns out to be a relevant operator which 
has the dimension smaller than 2. In order to see in which mode the relevant interaction 
opens the mass gap, we rewrite the interaction in abelian boson representation. From the 
formulas, J^r(l) ^ ^+(-)0c, J^r(l) 00 e^ 1 *^^ --^) , etc., and 



exp(— iv27T0c) exp(— iv27r < ) 
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etc., (5) 



where C = Cj l + 4>c,r and (p c = <fi Ct L — 4> c ,r Hl2| , the second term of eq.([|) becomes 
16\k cos v / 2tt(0c — 0/) cos ^/2n(p^\ Hence, the above relevant interaction may open ex- 
citation gaps in the spin mode ip = (0 C — 0/)/\/2 as well as in the charge mode <fyg). After 
integrating out the massive charge mode, <fy£\ we have the massive ip mode and the massless 
X = (0c + 0/)/V2 mode in the spin sector. This can be easily checked for v c = Vf where 
the 4> mode and the \ mode are decoupled. In the case with v c > Vf, by taking into account 
that Vf grows large toward the value Vf = v c by the renormalization, one can expect that 
the x mode is decoupled at the fixed point. Then, what happens for this massless spin mode 
X if we take into account the first term of eq.(§)? Notice that this interaction is marginal 
with the scaling dimension equal to 2. So, it is crucial whether it is marginally relevant or 
irrelevant. To see this, let us separate the interaction into two parts, 

H Zt + H int = ^k(Jc,L ■ Jf,L + Jc,R ■ J f,R) + A^(J C ,L " J f,R + J c,R ' J (6) 

The scaling equations for the couplings and \ r K can be obtained by expanding the 
partition function Z = tre - ^ in the Kondo coupling and by using the operator product 



expansion in SU(2) current algebra [|13| . [i4 



Up to the second order, the scaling equations read 

dAg dAfc (Ag) 2 

dlnL ' dlnL 2ir(v c + v f )' { ' 

We can see from (H) that the effect of H™ t is just to renormalize velocities, whereas H\ nt is 
marginally relevant, causing a spin gap. An important consequence of the marginal interac- 
tion is that the gap is given by the form of A s ~ exp(— bj \k) (b > 0), which is characterized 
by a Kosterlitz-Thouless-type transition. Therefore both of the charge and spin gaps open at 
half-filling, describing the fixed point of the Kondo insulator. These results agree with those 
of the numerical diagonalization and quantum Monte Carlo studies In particular, one 

can see that the spin gap of A s ~ exp(— b/Xx) deduced by numerical results || reflect the 
marginally relevant interaction in (§) 



We now turn to the case away from half-filling. Dropping irrelevant terms including 
oscillating factors in eq.([|), we end up with the interaction, 

Hm = X K J dx(J C:L (x) + J c ,r( x )) ■ S A x )i ( 9 ) 

which is still marginally relevant for the spin sector. In this case, charge excitations become 
massless, leading to a metallic state with the renormalized charge velocity. In order to see 
properties of the spin part, we first derive the scaling equation for the Kondo coupling Xk 
which is obtained similarly as in the half-filled case, 



d\ K (Xk) 2 



(10) 



dlnL 2tiv c 

We note that Caron and Bourbonnais also derived the scaling equations (|10D by using 
the Kadanoff- Wilson renormalization group method ||. In the infrared limit the Kondo 
coupling grows large towards the strong-coupling fixed point. We recall here that although 
in the — ► oo limit the ferromagnetic state appears ||, a paramagnetic state may be 
stabilized for a large parameter space in the metallic phase [15|]. Based on this observation, 



we assume the finite value of A^- for strong-coupling fixed point in the normal metallic phase. 
A simple physical picture at this fixed point is that the marginally relevant interaction (H) 
couples conduction electrons strongly with /-electrons to make singlet clouds, and remaining 
unpaired /-electrons form a massless spin mode which carries SU(2) currents. Strictly 
speaking, conduction electrons also contribute to the massless spin mode because of the 
finite value of A^ at the fixed point (see also discussions for the Fermi surface below). 
Therefore in the infrared limit, we have the effective hamiltonian 

H* = —J dx[{J c>L (x)} 2 + {J c , fi (x)} 2 ] + J dx[J ftL (x) ■ f f>L (x) + J ftR (x) ■ J f>R (x)}, 

(11) 

where v^ c ' and v j are the renormalized velocities of the charge and spin excitations. Note that 
the new SU(2) spin current j'j L ^ include both of conduction electrons and /-electrons. The 
above hamiltonian consists of the massless charge mode (holon) described by U(l) gaussian 



CFT and the massless spin mode (spinon) described by level-1 SU(2) CFT |L4]. Hence 
the fixed point in the metallic phase belongs to the universality class of the Tomonaga- 
Luttinger (TL) liquid as for the Hubbard model P,[16|-p0]l . Since this result is consistent 
with the conclusion deduced by Ueda et al. for A/ = |J, we believe that the introduction 
of the interaction A/ in eq.(l) may not change the essential physics of the Kondo lattice 
model even in the case away from half-filling. 

Now we wish to ask how is the volume of the Fermi surface? In ref . || , it has been claimed 
to be large. In the bosonization language, the Fermi surface may be described as follows. 
Since there is no relevant interaction in the charge sector, the left-going and right-going 
currents are decoupled and backward scatterings become irrelevant in the infrared limit. 
Thus the effect of the Kondo interaction on the charge sector at the strong coupling fixed 



point is just to give rise the tt/2 phase shift [21]. The total phase shift of conduction electrons 



due to Nf localized spins is nNf/2, which changes the " pseudo- Fermi surface" of holons 



from 2kp to 2kp + tt \\22\\ . This argument is not sufficient to deduce the real Fermi surface. 
As for the spin sector, the marginally relevant interaction @ again plays a crucial role 
to hybridize two kinds of spinons consisting of conduction electrons (pseudo-Fermi surface 
k-p) and /-electrons (vr/2), making the new " pseudo- Fermi surface" of massless spinons as 
kp + tt/2. Note that these pseudo- Fermi surfaces of holons and spinons are consistent with 
the formation of the charge and spin gaps at half-filling kp = tt/2. Combining these pseudo- 
Fermi surfaces, we can say that not only conduction electrons but also localized electrons 
contribute to the Fermi surface, forming the large Fermi surface, kp + tt/2, as we shall see 
more directly in the arguments for the single-particle Green function mentioned below. The 
result is in accordance with the conclusion of Ueda et al. ||. 

Let us now briefly discuss the critical exponents of correlation functions in the normal 
metallic phase. According to the above discussions, we have two kinds of primary fields 
<3>^ ± (x) (ho Ion) and ±(x) (spinon) with conformal dimensions |]l8l , |l9f , 



^ = \(i s -^±d})\ (12) 

which satisfy the requirements of U(l) and SU(2) current-algebra symmetries respectively. 
Here I c (J s ) and D c (D s ) are quantum numbers for the charge (spin) degrees of freedom, 
which obey the selection rules of Fermi statistics, D c = {I c + I s )/2 and D s = I c /2 mod 
1 [JTSIO. K p is a parameter which features U(l) critical line, depending on the bare pa- 
rameters Xx/t and A//i. The operators of conduction- and /-electrons can be expanded in 

terms of these primary fields as c a {x) = ^a±a± a A±A± < ^ > A ± ( a ')^ > A ± ^ a ')' e ^ c ' anc ^ hence vari- 
ous correlation functions in the asymptotic region are determined by those for the primary 
fields |18|]l9[| . For example, the single particle Green function for the conduction electrons, 
(c^(x)c cr (0)), is obtained by taking the quantum numbers (J c , I s , D c , D s ) = (1, 1, 0, 1/2). The 
resulting phase factor exp(2iD c khx) exp(2iD s k s x), where kh and k s are the "pseudo-Fermi 
surfaces" for holon and spinon given above, determines the position of the singularity in the 
Green function, i.e., the Fermi surface. Thus (cj.(x)c cr (0)) shows the singularity at k F + n/2 
(large Fermi surface) with the exponent rj = (K p + 1) 2 /AK P . Next we consider spin-spin 
correlation functions, (S£ (x, i)S~(0, 0)}, (Sf (x, t)Sj(0, 0)), and (S^(x, t)Sj (0, 0)), which 
show power-law decay with the same exponents determined by the primary field with 
(J c , I s , D c , D s ) = (0, 0, 1, -1/2) such that 

(St(x, 0)SJ (0, 0)) ~ ^■ + —(-l) x cos2k F x, etc., (13) 

where a = 1 + K p . We can deduce the correlation exponents in some limiting cases. For 
the weak coupling limit (A^ — > 0) the critical exponent may be characterized by that of free 
electrons with K p ~ 1. Near half-filling, the charge sector scales to the strong correlation 
limit characterized by the spinless-fermion exponent K p = 1/2. It is hence predicted that 
the spin exponent near half-filling takes a = 3/2 and that for the momentum distribution is 
9 = 1/8, for any value of the Kondo coupling A^. In general the correlation exponents may 



range 1/2 < K p < 1, similarly to the case of the Hubbard model [17|-|2"0|. It is instructive 



to note that the /-spin correlation (O) shows a power-law with 2k F + tc oscillation (neither 



2/cp nor tt). This long-distance behavior may be observed in the energy scale smaller than 
~ exp(— b/Xx), although the RKKY interaction may control the short- distance behavior, 



making characteristic structures around q = 2kp in correlation functions ||23|| . In particular, 
for small Xk the q = 2fcp structure may become more prominent [[4|,p4||. Even in this case, 
however, the low-energy critical behavior may be determined by the 2k-p + tt singularity so 
long as Xk is finite. 

Finally some comments are in order for bulk quantities. For example, one can expect the 
spin susceptibility to show a logarithmic low-temperature dependence in the metallic phase 
due to the leading irrelevant (marginal) operator such as Jf L ■ J t R , as observed in s = 1/2 



Heisenberg chain p5| , while it may exhibit exponential dependence at half-filling. Another 
remarkable feature is that as the electron density approaches half-filling at zero temperature, 
the charge susceptibility as well as the spin susceptibility may exhibit a divergence property 
reflecting the formation of the g observed for the charge susceptibility of the ID 

Hubbard model PS|| . Particularly, the spin susceptibility is predicted to behave like Xs — 
exp(— 6/A^)(l — n)^ 1 for small A^ reflecting the marginally relevant interaction (n: density 
of conduction electrons). 

In this paper, we have not discussed the effects of Xf in detail, which may appear through 
the renormalization of the velocity Vf. The scaling equation for Xf has a similar form to 
eq. (|TUD . So, there is a possibility that a massless spin phase may be stabilized at half- filling 
due to the large renormalization effect off/, if Xf exceeds a certain critical value A c initially. 
We think that the results obtained here, such as the large Fermi surface, are valid in the 
region for small Xf (< A c ). 

We have mainly discussed the case with the antiferromagnetic Kondo coupling. In the 
ferromagnetic Kondo case, it is known that a gap (a la Haldane gap) opens [4]. Unfortunately 
our formalism cannot describe the gap formation for this case. We think that another 
relevant mechanism which may be dropped in our formalism should be taken into account 
for generating the gap for the ferromagnetic case. This point should be clarified in the future 



S 



study. 
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